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PREFACE 

The  primary  purpose  of  this  report  is  to  fulfill  the  contract  require¬ 
ment  for  a  set  of  design  data  sheets,  similar  in  form  and  application 
to  Buships  Design  Data  Sheet  DDS  3801-2,  but  for  calculation  of  pres¬ 
sure  losses  at  high  gas  velocities  in  gas  turbine  ducting.  This  purpose 
is  accomplished  in  the  design  charts  and  design  procedures  found  in 
Sections  4  and  5.  These  have  been  prepared  in  such  form  that  they 
may  be  extracted  from  the  other  sections  of  the  report  to  form  a  self- 
contained  set  of  procedures  and  data  necessary  to  the  calculation  of 
high  speed  ducting  losses,  and  so  contain  the  necessary  elements  for 
a  Design  Data  Sheet. 

As  we  assembled  the  material  necessary  for  the  required  procedural 
and  data  sections,  we  realized  that  this  report  would  be  of  more  value 
to  those  who  are  ultimately  responsible  for  preparation  and  the  peri¬ 
odic  revision  and  improvement  of  the  Design  Data  Sheets  and  perhaps 
more  valuable  to  the  informed  user,  if  two  more  sections  were  to  be 

included  in  the  report: 

«** 

(1)  A  brief  introduction <o  the  parameters  and  general  methods  of  ~ 
compressible  (high  speed)  flow.  The  calculation  procedures, 
and  the  governing  flow  parameters,  are  in  most  cases  quite 
different  from  those  met -in  calculating  incompressible. (low 
speed)  duct  flows.  We  hope  that  the  short  introduction  included 
as  Section  2  of  the  report  will  be  helpful  to  those  who  are  primar¬ 
ily  familiar  with  the  methods  of  incompressible  flow  calculations, 
in  understanding  the  procedures  recommended  in  this  report,  and 
in  appreciating  some  of  the  reasons  for  their  comptaKlty  when 
compared  to  incompressible  methods. 
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(2)  Some  notes  on  the  data  sources  and  methods  that  were  employed 
In  the  preparation  of  the  design  charts  and  the  procedures.  Our 
survey  of  the  data  available  in  the  literature  on  compressibility 
(high  speed)  effects  in  duct  flows  showed  that  there  are  many 
areas  where  data  is  incomplete  or  inconclusive,  and  in  fact 
sometimes  conflicting.  This  is  also  true  to  a  lesser  extent  of 
the  available  data  on  incompressible  duct  flows  from  ’^us 
sources.  This  situation  resulted  in  the  necessity,  in  pre~>  ring 
the  final  charts  and  procedures  which  are  our  recommendations 
for  calculation,  of  applying  our  judgement  in  choosing  data  sources; 
of  generalizing  from  incomplete  data;  and  of  calculating  from  theo¬ 
retical  analysis,  to  extend  existing  data  or  methods  (e.  g. ,  the 
analysis  for  abrupt  area  decrease  in  Section  3.2.2). 

We  have  attempted  to  make  the  design  charts  and  calculation  proced¬ 
ures  as  complete  as  possible.  This  means  that,  in  some  cases,  where 
data  is  incomplete  or  inconclusive,  we  have  made  for  the  designer  rhe 
decisions  and  judgements  that  would  confront  him,  if  lie  were  working 
from  the  basic  data  available.  Under  these  circumstances,  we  have 
felt  it  would  be  useful  to  include  in  our  report  some  specific  indica- 
tions  of  the  sources  and  methods  which  led  to  our  recommendations. 

We  hope  that  these  notes  will  be  useful  if  a  critical  review  is  made  of 
these  methods  at  some  later  time  when,  presumably,  more  high  speed 
duct  flow  data  will  be  available  in  the  literature.  These  notes  on  methods 
and  sources  we  have  included  as  Section  3  of  the  report. 


Approved  by; 


James  R.  Turner 
Project  Manager 


L.  C.  Hoaglund 
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1.0  INTRODUCTION 


Shipboard  gas  turbine  inlet  and  exhaust  ducting  is  presently  designed 
for  low  speed  gas  flows,  of  the  order  of  100  ft/sec  and  less.  For 
these  velocities,  ducting  pressure  losses  may  be  estimated  with 
negligible  error  by  employing  the  assumption  that  the  gas  is  an  in¬ 
compressible  fluid.  This  assumption  implies  that  density  variations 
are  small  and  can  be  ignored  in  calculation  methods.  Pressure  loss 
calculation  methods  for  low  speed  flows  based  in  incompressible  as¬ 
sumptions  have  been  collected  and  formalized  in  Reference  1.  - 

If  it  is  desirable  to  design  inlet  and  exhaust  ducting  for  higher  flow 
velocities,  and  thus  obtain  ducting  of  smaller  diameter,  the  losses 
may  no  longer  be  estimated  by  the  assumptions  of  incompressible 
flow.  At  high  subsonic  gas  velocities,  above  about  Mach  Number  0.  2, 
the  magnitude  of  the  density  changes  accompanying  the  pressure 
changes  due  to  accelerations,  decelerations,  and  losses  should  not  be 
ignored  in  calculation.  Furthermore,  the  loss  coefficient  character¬ 
izing  the  pressure  loss  in  each  ducting  component,  usually  expressed 
as  a  fixed  percentage  of  the  velocity  head  for  a  component  of  given 
configuration  in  incompressible  flow,  generally  depends  also  on  the 
flow  velocity  (or  Mach  Number)  in  high  speed  flows.  The  presenta¬ 
tion  of  a  calculation  procedure  for  high  speed  duct  flows  therefore 
requires  two  extensions  of  low-  speed  calculation  procedures: 


1.  Formulation  of  a  calculation  procedure  which  will  take  into 
account  the  effects  of  compressibility  on  local  values  of  flow 
parameters  (pressure,  density,  velocity,  etc.) 


2.  Correlation  of  available  data  or.  the  «j er  of  oar./.iv 
on  the  pressure  losses,  or  loss  coefficients  which  cha 
the  pressure  losses,  in  various  due:  components. 


This  report  describes  a  calculation  system  which  permits  the  cal¬ 
culation  of  ducting  losses  for  high  speed  flows.  In  the  following 
sections  the  methods  and  data  sources  employed  are  described, 
and  a  system  is  devised  for  application  to  duct  pressure  loss  cal¬ 
culations.  This  system  has  been  reduced  to  a  series  of  charts  for 
computational  aid.  These  charts  are  included  in  a  separate  section 
of  this  report  and  are  arranged  so  that  they  may  be  extracted  from 
the  report  to  form  a  self-contained  calculation  system  similar  in 
form  and  intent  to  the  material  of  Reference  1.  Although  the  sys¬ 
tem  is  not  as  simple  in  form  and  application  as  the  incompressible 
procedures,  they  represent,  in  our  opinion,  the  simplest  calculation 
procedures  which  can  be  applied  to  high  speed  flow  calculations. 
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2. 0  INTRODUCTION  TO  PARAMETERS  AND  GENERAL  METHODS 
OF  COMPRESSIONS  FLOW 

Accounting  for  gas  compressibility  effects  in  the  calculation  of  ducting 
pressure  losses,  necessary  for  high  speed  duct  flows,  introduces 
several  new  variables  into  the  calculations.  Rather  extensive  changes 
to  incompressible  flow  methods  and  procedures  must  be  made.  The 
purpose  of  this  section  is  to  provide  a  brief  outline  of  parameters  and 
methods  employed  in  this  report  which  are  peculiar  to  compressible 
flow  calculations.  Those  analytical  relations  developed  which  are 
important  to  the  calculation  procedure  have  been  reduced  to  chart 
form  in  later  sections. 

2. 1  Mach  Number 

The  Mach  Number  is  the  basic  measure  of  the  relative  importance 
of  compressibility  effects.  The  Mach  Number  is  the  ratio  of  the  fluid 
velocity  to  the  local  speed  of  sound: 


where 

V  =  fluid  velocity 

c  =  local  speed  of  sound  =  y gQ  k  RT 

Only  subsonic  duct  velocities  (i.  e. ,  M  <1. 0)  have  been  considered. 
For  Mach  Number  M  <0.2,  compressibility  effects  are  small,  and 
incompressible  calculation  methods  can  be  applied  with  sufficient 
accuracy. 


Shapiro  (Reference  2)  and  others  have  shown  that  most  compressible 
flow  calculations  are  most  easily  carried  out  if  the  calculation  meth¬ 
ods  are  formulated  in  such  a  way  that  Mach  Number  is  the  independent 
variable  of  the  calculation  procedure.  For  many  of  the  pressure-loss 
calculations  which  follow,  the  Mach  Number  has  therefore  been  chosen 
as  the  independent  variable.  Generally,  losses  in  ducting  components 
which  do  not  involve  area,  change  cause  the  Mach  Number  to  increase. 
Therefore,  it  is  possible  that  in  a.  long  ducting  system,  incompressi¬ 
ble  calculation  methods  may  be  appropriate  in  sections  near  the  inlet, 
with  compressible  methods  becoming  necessary  as  the  Mach  Number 
increases. 


2.  2  Total  Pressure  and  Loss  Coefficients. 


The  total,  or  stagnation,  pressure  is  defined  as  the  pressure  which 
a  gas  stream  would  reach  if  it  were  decelerated  to  zero  velocity 
without  incurring  any  losses  in  the  deceleration.  The  total  pressure 
is  the  pressure  measured  by  a  Pitot  tube,  or  an  "impact"  probe.  For 
incompressible  flow,  the  total  pressure  is  the  sum  of  the  static  pres¬ 
sure  and  the  velocity  head. 


=  P  + 


P  V' 


So 


(2.1) 


where 

pQ  =  total  pressure  (impact  probe  pressure) 

p  =  static  pressure  (pressure  at  a  wall  pressure  tap) 

p  =  fluid  density 

V  =  fluid  velocity 
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The  change  in  total  pressure  in  a  duct  flow  is  a  measure  of  the  losses 
in  the  flow.  The  losses  may  be  characterized  by  a  dimensionless  loss 
coefficient,  Cpo> 


cp0  ■ 


;  oi 


1 

2 


(2-2) 


where  1  and  2  represent  an  upstream  and  downstream  station,  re¬ 
spectively.  For  incompressible  flow  in  a  constant  area  duct, 


>iv:  . 

2  So  2  So 

and  so, 

Pol  “  Po2  =  P1  “  p2 


(2.3) 


(2.4) 


That  is,  the  difference  in  total  pressure  between  two  stations  in  the 
duct  is  equal  to  the  difference  in  static  pressure,  so  that  the  flow 
losses  are  described  equally  well  by  either  total  pressure  change  or 
static  pressure  change.  'This,  however,  is  not  true  for  a  compres¬ 
sible  (i.  e. ,  high  speed)  flow. 


For  a  compressible  flow,  the  simple  relation  between  total  pressure, 
static  pressure,  density,  and  stream  velocity  given  by  Equation  2. 1 
cannot  be  applied  directly.  This  equation  does  not  apply  because  the 
fluid  density  will  change  during  the  deceleration;  that  is,  the  density, 
p,  in  Equation  2. 1  is  not  a  constant.  The  relation  corresponding  to 
'-Equation  2. 1  for  compressible  flow  is  best  expressed  in  terms  of  die 
Mach  Number: 
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(1 


k  -  1 
+  ~~2 — 


9  k  -  i/k 
M  ) 


(2.5) 


where  k  =  ratio  of  specific  heats,  a  gas  property. 

The  losses  can  be  expressed  as  a  change  in  total  pressure  in  compres¬ 
sible  flow;  this  change  is  not  equivalent  to  the  static  pressure  change, 
as  in  incompressible  flow.  The  definition  of  loss  coefficient,  Equation 
2.2,  may  be  rearranged  to  more  convenient  form  through  the  relation 


so  that. 


Pox  "  P02 


TkPl  Mj2 


Combined  with  Equation  2.  5  this  becomes. 


(2.6) 


(2.7) 


1  -  P°2 
Pox 


k  -  i 
~T~ 


(2-  8) 


Thus,  if  the  loss  coefficient,  CpQ,  and  the  inlet  conditions,  and  pQi , 
are  known  for  a  given  ducting  component,  the  exit  total  pressure  for 
that  component  can  be  determined  from  equation  2.8;  the  exit  Mach 
number,  M2,  is  determined  from  continuity  relations,  as  will  be  shown 
in  Section  2.4.  These  exit  conditions  will,  of  course.  .. •  .«ic  1/. 
conditions  for  the  following  uucting  component. 
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2.3  Total  Temperature 

The  total,  or  stagnation,  temperature  is  similar  to  the  total  pressure 
in  definition.  That  is,  it  is  the  temperature  which  a  gas  stream  would 
reach  if  it  were  decelerated  to  zero  velocity.  It  is  best  expressed  in 
terms  of  the  Mach  Number: 

—  =  i  +  111  M2  (2.9) 

T  2 

where 

T  =  total  temperature  (temperature  of  stream 
decelerated  to  zero  velocity) 

T  =  static  temperature  (temperature  of  the 
flowing  gas  stream) 

Unlike  the  total  pressure,  however,  the  total  temperature  is  un- 
'affected  by  losses.  In  a  duct  flow,  the  total  temperature  can  change 
only  through  heat  transfer  to  or  from  the  gas  stream.  Note,  however, 
that  the  static  temperature  may  vary  considerably,  as  by  acceleration 
or  deceleration  of  the  gab  stream,  according  to  Equation  2.9. 

If  we  assume  that  heat  transfer  to  the  surroundings  is  negligibly  small 
in  a  duct  flow  (i.e. ,  the  flow  is  adiabatic),  the  total  temperature  re¬ 
mains  constant,  irrespective  of  the  flow  losses.  The  assumption  of 
adiabatic  flow  has  been  made  in  the  calculation  methods  presented  in 
this  report.  Therefore,  for  a  gas  turbine  intake  duct  flow,  the  total 
temperature  is  the  temperature  of  the  atmosphere  (i.e.,  the  tempera¬ 
ture  of  the  inlet  air  when  it  was  at  zero  velocity),  and  is  constant 
throughout  the  inlet  ducting;  for  a  gas  turbine  exhaust  duct  flew,  the 
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total  temperature  is  determined  by  the  exhaust  condition  at  the  turbine, 
and  is  constant  throughout  the  exhaust  ducting. 

The  assumption  of  adiabatic  flow  in  the  intake  duct  is  quite  accurate, 
since  only  small  temperature  differences  exist  between  the  duct  flow 
and  the  surroundings.  For  the  hot  exhaust  flow,  the  adiabatic  condi¬ 
tion  is  not  so  nearly  achieved,  though  insignificant  error  is  probably 
introduced  for  insulated  ducting.  To  account  for  the  effects  of  heat 
transfer  to  the  surroundings  in  the  exhaust  duct  flow  would  preclude 
any  simple  calculation  method  for  the  exhaust  duct  flow  losses. 

2.4  Mass  Flow  Parameter  and  the  Continuity  Relation 

For  incompressible  flow,  the  continuity  relation  may  be  expressed  as 

» 

m  =  .pA^V^  e  p  A ^2  (2.10) 

Since  density  is  constant  between  any  two  sections,  the  velocity  will 
change  only  by  virtue  of  nr.  area  change.  In  a  compressible. flow, 
however,  velocity  changes  (and  thus  changes  in  Mach  Number)  car.  ' 
occur  in  a  constant  area  due:,  because  of  changes  in  density  resulting 
from  pressure  losses,  ”hu.s,  the  continuity  relation  must  be  applied 
along  with  the  loss  relations  to  determine  Mach  Number  changes  ac¬ 
companying  ducting  component  losses.  For  compressible  flow,  the 
continuity  relation  may  be  conveniently  expressed  as 

=  0. 532  F  (M)  (2.11) 

P  A 
TO 

where  the  constant  0.  532  applies  for  a  gas  of  k  =  1.  4  and  molecular 
weight  of  29.  The  function  of  Mach  Number,  F  (M),  is 
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This  function  is  the  reciprocal  of  the  function  A/A*  tabulated  in  Table 
30  of  Reference  3. 

For  an  adiabatic  duct  flow  (T  =  constant),  the  conBnuIty~equation  be¬ 
comes,  with  =  m^, 

P0l''i  F(Mt)  =  po2  A2F<M2)  (2.13) 

Thus,,  for  a  given  duct  component,  with  fixed  area  ratio,  the  exit 
Mach  Number,  M2,  may  be  found  if  inlet  conditions  pQi  and  are 
known,  and  the  ratio  P02/P0i  is  known  from  loss  coefficient  data 
(see  Section  2.2). 

Equation  2.11  gives  a  convenient  way  to  determine  the  Mach  Number 
for  a  duct  flow  when  mass  flow,  area,  total  pressure,  and  total  tem¬ 
perature  are  known. 


The  term 


is  denoted  the  "mass  flow  parameter”;  for  a  given  gas  this  parameter 
depends  on  the  Mach  Number  only.  A  similar  parameter  employing 
static  pressure, 


P  A 


can  also  be  determined  as  a  function  of  Mach  Number.  These  param¬ 
eters  are  plotted  vs.  Mach  Number  in  a  later  section,  to  give  a  simple 
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computational  aid  in  determining  Mach  Number  when  other  flow  con¬ 
ditions  are  specified. 

5  Choking  in  Compressible  Flow 

The  mass  flow  rate  in  a  duct  flow  of  compressible  gas  is  limited  by 
the  phenomenon  of  ’'choking'1.  This  may  be  explained  by  the  following 
considerations. 


Consider  a  simple  ducting  component  which  has  a  continuous  area 
change  (e.g. ,  a  nozzle),  with  specified  areas 


and  and  flow  at  to.al  temperature  T  .  If  we  assume  that  there 
are  no  losses,  pQi  =  pq2  *P0-  Let  us  first  hold  the  upstream  total 
pressure,  pQ^,  constant  and  continuously  reduce  the  downstream  static 
pressure,  p2<  From  Equation  2.i 


fo2 
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.  0  k-l/k 

m2  ) 


(2.14) 


the  Mach  Number  M2  will  increase  as  static  pressure  p2  is  reduced. 
The  effect  on  mass  flow  rate  is  determined  through  Equation  2. 11, 


=  0.532  F(M2) 


(2. 15) 


and  Equation  2. 12, 


PoA2 


F(M) 


M 


k  +  l/2(k-l) 


(2.16) 


Sketch  1  shows  the  variation  of  F  (M)  with  Mach  Number.  As  M  is 
increased,  F  (M)  increases  to 


a  maximum  value  at  M  =.1.0,  and  then  decreases.  Thus,  for  our  ex¬ 
ample,  as  the  pressure  is  decreased,  M2  will  increase  and  the  mass 
flow  rate  will  increase  until  a  Mach  Number  M2  =1.0  is  reached,  at 
which  condition  the  flow  wili  be  a  maximum.  This  occurs  at  a  value  of 
8g^/p^  =  0.  528,  for  a  gas  with  k  =  1.  4.  It  can  be  shown  that  the  Mach 
0"JL  Number  cannot  pass  through  M2  =1.0  to  supersonic  values  in  the  min¬ 
imum  area  A2*  .  Thus,  as  the  pressure  p2  is  further  reduced,  the 
value  of  M2  remains  at  unity,  and  the  mass  flow  rate  remains  constant. 
This  maximum  value  of  flow  rate  is  the  "choked"  flow  rate.  Thy 


♦See  Reference  2,  Chapter  4  for  a  rigorous  explanation  of  this  limitation. 
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variation  of  flow  rate  with  downstream  static  pressure,  pj,  at  a 
constant  upstream  total  pressure  pQ  is  thus  shown  by  Sketch  2. 


Note  that  when  the  choked  flow  condition  (M2  *  1. 0)  is  attained,  the 
inlet  Mach  Number,  M^,  also  has  a  definite  value  which  is  fixed  by 
the  area  ratio  From  Equation  2. 13,  at  M2  c  1.0,  F(M2  *1.0) 

«1.0  and  Po2'Pol'P0 


F(MX) 
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Therefore, 

N 


F<M1> 


(2.17) 


(2. 18) 


Thus,  for  this  duct  component,  there  is  a  limiting  value  of  inlet  Mach 
Number,  M^,  which  cannot  be  exceeded  due  to  the  "choking"  phenome¬ 
non. 


If  we  now  assume  that  there  is  a  total  pressure  loss  in  the  above  nozzle, 
then  P02/P0i  <  1-0,  and  the  inlet  Mach  Number  at  the  choked  condition 


is  given  by 


F  (Ma  ) 


po2  A2 
pol  A1 


(2.19) 


Reference  to  Sketch  1  will  show  that  the  limiting  inlet  Mach  Number 
will  now  be  lower  than  for  the  nozzle  without  losses,  since  F  (Mj, ) 
is  smaller  by  the  ratio  P^/P^*  Thus,  the  limiting  inlet  Mach 
Number  due  to  choking  for  a  duct  component  is  influenced  (1)  by 
area  change  and  (2)  by  losses;  either  of  these  influences  may  cause 
a  choked  condition  at  the  exit,  with  an  inlet  Mach  Number  less  than 
unity.  For  example,  a  straight  constant- area  duct  with  frictional 
losses  will  choke  at  an  inlet  Mach  Number  less  than  unity,  with  the 
choking  value  of  inlet  Mach  Number  depending  on  the  duct  length, 
diameter,  and  wall  roughness;  a  constant- area  bend  likewise  will 
choke  at  an  inlet  Mach  Number  less  than  unity  because  of  losses  in 
the  bend. 

The  practical  significance  ;o  :.-.e  duct  system  designer-can.be  illus¬ 
trated  by  an  example.  Suppose  that  we  wish  to  design  a  minimum 
size,  constant  area  inlet  rir  duct  for  a  gas  turbine.  The  air  is  to 
be  drawn  from  atmosphe:  w  conditions,  so  that  pQi  and  Tq  are  fixed; 
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the  required  turbine  mass  flow,  m,  will  be  specified.  Equation  2. 11 
gives 


mlTf“ 

- L_iL  =  0.532  F(M) 

PoA 


(2. 11) 


Let  us  suppose  first  that  the  inlet  duct  is  so  short  that  losses  will  be 


completely  negligible.  Then  pQ,  Tq,  and  m  are  all  specified,  and 
the  minimum  area  duct  will  occur  at  the  maximum  value  of  the  mass 


flow  parameter,  myT^/poA  .  This  corresponds  to  the  maximum 
value  of  F(M),  which  in  turn  occurs  at  M  =1  (also  F(M)  =  1.0),  the 


choked  condition.  Thus,  the  minimum  area  duct  is  determined  by 


0.532 


(2.20) 


when  the  specified  values  of  m  ,  T  ,  and  p  are  inserted.  If  a  duct 

o  o 

any  smaller  than  this  were  used,  the  specified  mass  flow  would  not 
be  attained,  regardless  of  how  low  the  gas  turbine  compressor  re¬ 
duced  the  pressure  at  the  duct  exit. 


If  we  now  allow  a  long  inlet  duct  in  which  losses  are  not  negligible, 

(and  F(M2)  cannot 
exceed  unity),  the  maximum  mass  flow  parameter  becomes 


then  po2  <  pQ^.  Since  M2  cannot  exceed  unity 


Since  pq2  is  less  than  pQi,  the  duct  area  required  to  admit  the  spec¬ 
ified  flow  becomes  larger  than  for  the  short  inlet  duct,  because  of 
the  duct  losses.  The  limiting  inlet  Mach  Number  will  1>L  amnllor  than 
unity. 
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Near  choking,  duct  component  losses  generally  become  so  large  that 
^■"^it  is  not  practical  to  design  for  choked  or  near-choked  ducting.  How¬ 
ever,  it  is  important  to  realize  that  regardless  of  how  much  pressure 
loss  is  allowable,  there  are  minimum  duct  areas  that  will  admit  a 
specified  mass  flow  at  fixed  inlet  conditions  (e.g. ,  atmospheric  con¬ 
ditions). 

2.6  Gas  Properties 

»• 

All  the  relations  of  the  calculation  methods  presented  assume  that 
the  perfect  gas  equation  of  state, 

P  =  P  R  T 

is  valid  for  the  duct  flows  considered.  In  addition,  the  charts  have 
been  prepared  for  a  gas  with  ratio  of  specific  heats  of  1.4,  and  mo¬ 
lecular  weight  29.  These  values  will  apply  with  sufficient  accuracy 
to  either  atmospheric  inlet  air  or  gas  turbine  exhaust  gases. 


9 
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3.  0  SOURCES  AND  METHODS  EMPLOYED  IN  DEVELOPMENT  OF 
CALCULATION  PROCEDURES  FOR  DUCT  LOSS  CALCULATION 

The  data  sources  and  correlation  methods  on  which  the  calculation 
procedures  presented  in  this  report  are  based  are  outlined  in  this 
section.  In  some  cases  the  available  data  is  incomplete,  so  that 
approximations  and  extrapolations  have  been  necessary;  in  some 
cases,  data  from  various  sources  conflict,  so  that  some  judgment  - 
has  been  necessary  in  choosing  data  for  procedures.  Under  these 
circumstances,  it  seems  desirable  to  supplement  the  recommended 
calculation  procedures  with  this  outline  of  sources  and  methods. 

3. 1  Straight  Constant  Area  Ducts 

The  procedures  for  straight  ducts  are  based  on  the  calculation  methods 
developed  by  Shapiro,  and  others.  These  are  developed  in  detail  in 
Reference  2,  and  are  a  standard  method  for  compressible  flow  calcu¬ 
lation  to  include  the  effects  of  friction.  The  necessary  compressible 
flow  parameters  for  this  calculation  are  plotted  in  Figure  D-7. 

Duct  friction  factors  are  taken  from  Moody,  Reference  4,  and  charts 
are  reproduced  from  thiS  reference  as  Figures  D-3  and  D-4.  Keenan 
and  Neumann,  Reference  5,  show  that  there  is  no  significant  effect 
of  compressibility  on  friction  factor  for  subsonic  flows,  so  that  the 
friction  factors  of  Moody  may  be  applied  without  Mach  Number  correc¬ 
tion. 

The  friction  factors  of  Moody  were  determined  in  fully- developed  duct 
flow's  and  do  not  apply  to  flows  in  inlet  lengths  of  ducting,  within  about 
30  to  50  diameters  of  the  duct  entrance.  In  this  entrance  length,  ve¬ 
locity  profiles  are  not  yet  fully  developed.  Reference  7  shows  that,  for 
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an  inlet  length  with  a  smooth  entrance  (e.g. ,  a  bellmouth),  the  flew 
is  laminar  until  a  length  Reynolds  Number,  Re  r  «  10^,  is  attained 
(Re^  =p  V x/jx ,  x  being  the  distance  from  the  duct  entrance).  The 
theory  of  Langhaar,  Reference  6,  has  been  shown  by  Shapiro  to  pre¬ 
dict  the  correct  value  of  friction  factor  in  this  laminar  region.  Curves 
prepared  from  Langhaar' s  theory,  taken  from  Reference  7,  are  re¬ 
produced  as  Figure  D-5.  At  the  duct  length  corresponding  to  Re  =10^, 
a  transition  to  turbulence  occurs;  Reference  7  presented  data  for 
friction  factors  in  this  turbulent  region.  The  curves  of  Figure  D-6 
represent  an  average  friction  factor  (calculated  from  curves  faired 
through  this  data)  to  be  applied  in  this  turbulent  inlet  length.  The 
friction  factor  in  the  turbulent  inlet  region  is  shown  as  a  ratio  to  the 

friction  factor  for  fullv-develoced  turbulent  flow  at  the  same  duct 

¥  «. 

Reynolds  Number  (Re  *=pVD/(i);  this  ratio  is  larger  than  unity. 
Therefore,  losses  are  larger  in  an  inlet  length  than  in  a  fully-developed 
flow.  If  the  duct  entrance  is  not  smooth,  or  if  there  is  a  screen  in  the 
entrance,  Shapiro  has  shown  that  there  is  no  laminar  inlet  length,  but 
that  the  entire  inlet  flow  length  is  turbulent.  Thus,  for  such  entrance 
conditions,  the  tur-bulent  inlet  friction  factors  of  Figure  D-6  should  be 
applied  to  the  entire  inlei  length. 

Calculation  procedures  have  been  devised  which  take  into  account 
these  special  cases  of  inlet  flow  lengths  and  determine  the  appropriate 
friction  factors  according  to  the  above  considerations.  These  are 
presented  as  procedures  P-9  and  P-iO.  Inspection  of  these  procedures 
will  show  that  they  are  complex  and  lengthy.  For  this  reason,  a  simpler 
approximate  procedure  for  inlet  length  calculations  has  been  prepared, 
procedure  P-3.  This  approximate  procedure  employs  the  ful’v 
flow  friction  factors  of  Moody,  and  ■'ill  u:  u -•rcstin ti.y  !>)•  ■[>,;, 
losses.  It  has  been  recommended  for  use  in  the  duct  calculation  method, 
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however,  because  the  improved  accuracy  of  the  more  complex  methods 
is  not  generally  warranted  for  duct  system  calculations,  because  cf 
the  uncertainties  in  the  loss  calculations  cf  other  duct  system  compon¬ 
ents.  If,  however,  a  duct  system  consists  primarily  of  straight  runs 
of  ducting,  or  if  the  designer  judges  the  improved  accuracy  desirable, 
the  more  accurate  procedures  P-9  and  P-10  can  be  used. 

The  losses  for  beilmcuth  entrances  incorporated  in  these  procedures 
have  been  estimated  from  the  data  cf  Henry,  Reference  8.  The  losses 
for  sharp- edged  inlets,  and  for  inlets  with  screens  and  gratings,  have 
been  computed  by  a  method  similar  to  that  described  in  the  section  on 
abrupt  area  decrease,  incorporating  an  extrapolation  of  the  data  for 
contraction  coefficients,  Reference  12. 

3. 2  Area  Change 

3. 2. 1  Abrupt  Area  Increase 

A  calculation  method  for  prediction  of  losses  in  subsonic  flows  through 
an  abrupt  area  increase  ("sudden  expansion")  has  been  developed  by 
Hall  and  Orme,  Reference  9.  These  authors  have  verified  their  cal- 
culatior.s  by  experiment  over  a  range  cf  area  ratios  from  4  to  30,  at 
entering  Mach  Numbers  from  0.25  to  i.0.  Cole  and  Mills,  Reference 
10,  have  confirmed  the  method  by  ex:  eriment  in  the  area  ratio  range 
from  i  to  4.  This  analytical  method  has.  boor,  applied  here  to  the  cal¬ 
culation  cf  total  pressure  losses  and  Mach  Number  change,  and  results 
are  plotted  in  the  design  charts  of  Figure  D-9. 
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The  experimental  data  of  Hall  and  Orme  indicate  that  the  total  pressure 
loss  and  Mach  Number  change  take  place  in  a  mixing  length  of  4  diame¬ 
ters  downstream  of  the  sudden  area  change.  In  the  calculation  proced¬ 
ure,  the  total  pressure  loss  attributed  to  the  aiea'  change  is  therefore 
taken  to  include  the  loss  of  4  diameters  of  downstream  duct.  In  the 
calculation  of  losses  in  other  ducting  components  in  which  the  loss  is 
attributable  to  a  separation  and  subsequent  mixing  of  the  flow  (e.  g. 
sudden  contractions,  screens,  sharp- edged  entrances,  etc.),  the  total 
pressure  loss  and  Mach  Number  change  attributed  to  that  component 
have  likewise  been  taken  to  include  4  downstream  duct  diameters,  on 
the  basis  of  the  experimental  evidence  of  Hall  and  Orme  (and  others, 
e.  g. ,  Reference  13),  unless  some  data  to  the  contrary  has  been  avail¬ 
able  for  that  component. 

3.2.2  Abrupt  Area  Decrease 

Apparently,  the  problem  of  compressible  flow  through  abrupt  area  - 
decreases  has  not  been  previously  considered  in  the  open  literature. 

An  analysis  has  therefore  been  developed  according  to  the  model  for 
incompressible  flow  suggested  in  Reference  11. 


/T 


Sketch  I 
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As  shown  in  Sketch  1,  the  flow  is  assumed  to  separate  from  the 
sharp  corner  of  the  area  change  to  produce  a  free  streamline  flow, 
contracting  to  the  streamtube  area,  A£.  Turbulent  mixing  then 
occurs  between  sections  2"  and  2  until  the  flow  fills  the  area  A2» 
at  essentially  constant  velocity  across  the  duct  cross-section.  The 
flow  from  section  2”  to  2  is  treated  as  an  abrupt  area  increase,  as 
described  in  Section  3. 2. 1.  The  area  Aj'  must  be  related  to  the 
known  areas,  A^  and  A2,  by  experiment,  since  there  is  no  analy¬ 
tical  method  for  prediction.  Values  of  area  ratio  A^/A^  have  been 
taken  from  the  data  correlations  of  Reference  12,  which_gives  the 
area  ratio  Ag  / A ^  as  a  function  of  A2/A^  and  inlet  Mach  Number,  M^. 

The  results  of  this  analysis  are  presented  in  Figure  D-10  which 
shows  exit  Mach  Number  and  total  pressure  ratio  as  a  function  of 
area  ratio  and  irtlet  Mach  Number.  Total  pressure  loss  and  Mach 
Number  change  include  4  diameters  of  downstream  ducting.  Although 
pressure  losses  are  small,  note  that  moderate  area  changes  produce 
large  increases  in  Mach  Number,  so  that  only  small  area  decreases 
can  be  afforded  in  high  speed  flows. 

3.  2.  3  Diffusers 

Although  a  large  amount  of  literature  describing  diffuser  investiga¬ 
tions  is  available,  correlations  of  these  data  have  as  yet  produced 
no  systems  for  diffuser  design  or  loss  calculation  that  are  entirely 
satisfactory.  Available  correlations  of  performance  have  been 
formulated  on  the  basis  of  diffuser  geometry  alone.  Various  inves¬ 
tigations  have  shown,  however,  that  there  can  be  substantial  effects 
of  Reynolds  Number,  inlet  boundary  layer  thickness  (i.  e. ,  inlet 
velocity  profile),  compressibility,  and  discharge  duct  geometry. 
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However,  there  is  not  sufficient  data  yet  available  to  include  such 
effects  in  a  general  correlation,  so  that  we  must  rely  on  correlations 
of  performance  with  geometry  alone.  Such  correlations  exhibit  con¬ 
siderable  scatter  of  the  data,  because  of  the  influence  of  the  effects 
noted  above,  and  so  are  subject  to  some  uncertainty  in  calculation. 

The  most  complete  correlation  available  in  the  open  literature  is  that 
of  Patterson,  Reference  13.  The  design  curves  of  design  chart  D- 11 
of  this  report,  taken  from  the  correlations  of  Patterson,  are  recom¬ 
mended  for  determination  of  incompressible  flow  losses  in  conical 
diffusers. 


Various  schemes  have  been  proposed  for  correlating  rectangular 
diffusers  with  conical;  the  method  recommended  here  is  to  assume 
the  loss  in  a  diffuser  of  rectangular  cross-section  to  be  that  of  an 


equivalent  conical  diffuser  circumscribed  about  the  rectangular  dif¬ 


fuser*  as  in  the  following  sketch. 
RECTANGULAR 


I 
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High  speed  flow  effects  on  diffuser  performance  are  shown  by  Little 
and  Wilbur,  Reference  15,  and  by  Young  and  Green,  Reference  14, 
for  several  conical  and  rectangular  diffusers.  Correlations  from 
the  data  of  these  investigators  have  been  made,  and  are  shown  on 
design  chart  D-ll  for  use  as  correction  factors  to  the  incompres¬ 
sible  loss  coefficients.  These  correlation  curves  are  based  on  data 
from  only  a  few  diffuser  geometries,  and  should  be  considered  to 
be  subject  to  revision  if  more  data  becomes  available  in  the  future. 
The  effects  of  increasing  Mach  Number  on  performance  of  a.  diffuser 
with  a  fixed  inlet  boundary  layer  condition  has  been  shown  in  Refer¬ 
ence  15  to  be  smaller  than  the  effect  of  increasing  inlet  boundary 
thickness  at  a  fixed  Mach  Number.  Such  results  as  these  indicate 
that  there  will  be  considerable  uncertainty  in  diffuser  loss  calcula¬ 
tions.  The  recommendations  in  the  curves  of  D-ll  simply  represent 
the  best  correlations  presently  possible  from  existing  data  sources. 
Referring  to  the  Mach  Number  correction  curves,  note  that  the  total 
pressure  loss  for  conical  diffusers  begins  to  increase  markedly  at 
entering  Mach  Numbers  of  0.  4  to  0.  5;  choking  occurs  in  the  range  of 
0.  6  to  0.7  Although  rectangular  d?"isers  do  not  show  as  marked  an 
^increase  in  loss  with  increasing  Mach  ..amber,  the  incompressible 
losses  are  much  larger  than  for  conical  diffusers;  thus,  the  total  loss 
always  is  larger  than  for  comparable  conical  diffusers. 

The  following  remarks  will  serve  as  a  guide  to  the  designer: 

1.  An  angle  between  diverging  walls  of  about  6°  - 10°  gives  best 
diffuser  performance. 

2.  For  a  given  area  change  per  unit  length  of  duct,  diffusers  at 
circular  cross-section  give  the  best  effectiveness,  with  square  cross- 
sections  next. 
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3.  The  inlet  velocity  distribution  has  been  shown  to  affect  diffuser 
performance  (e.  g.,  Reference  15).  A  velocity  distribution  with  a 
thin  boundary  layer  at  diffuser  entrance  gives  better  diffuser  per¬ 
formance  than  a  thick  entry  boundary  layer,  with  losjes  differing 
by  as  much  as  200%  between  the  two  cases.  Although  little  can  be 
done  in  a  ducting  system  to  control  boundary  layer  thickness,  this 
effect  indicates  that  better  diffuser  performance  can  be  expected 
for  a  diffuser  relatively  near  the  duct  system  inlet  (less  than  30 
diameters)  than  would  be  expected  after  sufficient  length  for  fully- 
developed  flow  to  be  established;  and  that  a  diffuser  following  a  bend, 
or  some  similar  component  which  may  generate  a  thick  or  separated 
boundary  layer,  will  have  relatively  poorer  performance. 

4.  Complete  pressure  recovery  is  not  realized  at  the  exit  of  the 
diffuser;  about  4  to  6  exit  duct  diameters  downstream  of  the  exit  are* 
required  for  complete  recovery. 

5.  For  most  commonly  used  area  ratios,  (up  to  about  6:1)  an  abrupt 
area  increase  results  in  less  total  pressure  loss  than  a  diffuser  with 
wall  angle,  2  $  greater  tKan  40°  -  50°. 

6.  Total  pressure  loss  increases  with  increasing  inlet  Mach  Number; 
the  effect  is  more  severe  with  large  wall  angles  and  also  with  thick 
entry  boundary  layers. 

7.  Various  schemes  such  as  boundary  layer  control  by  suction, 
injection,  and  vortex  traps,  and  delay  of  flow  separation  by  curved 
walls,  have  been  devised  to  improve  diffuser  performance.  Vanes 
have  been  used  to  stabilize  the  flow  and  thus  improve  performance. 


These  possibilities  have  not  been  considered  in  this  study,  since  no 
definitive  design  rules  can  be  laid  down  and  the  final  design  must  re¬ 
sult  from  testing  to  determine  empirically  an  optimum  configuir.uon. 

/ 

8.  When  space  limitations  would  require  a  diffuser  with  large  wall 
angle,  it  is  often  more  desirable  to  employ  an  efficient  diffuser  (wall 
angle,  2  g,  of  about  10°)  followed  by  an  abrupt  area  increase  to  the 
desired  area. 

3.2.4  Converging  Duct  (Nozzle) 

Losses  in  converging  ducting  are  small  compared  to  those  in  diverging 
ducting  (diffusers)  because  the  accelerating  flow  in  a  converging  duct 
presents  a  favorable  pressure  gradient  for  the  wall  boundary  layer. 

For  small  cone  angles  (i.  e. ,  less  than  30°  included  angle)  there  are 
generally  no  separ3tion-mixing  losses.  There  is  little  information 
available  in  the  literature  on  Mach  Number  effects  in  converging  sec¬ 
tions  with  downstream  ducting;  most  available  data  deals  with  nozzles 
exhausting  to  a  large  ch<  mber  or  atmosphere.  Therefore,  it  is  recom 
mended  that  converging  ducting  with  included  cone  angle  smaller  than 
30°  be  treated  as  straight  ducting,  with  L/D  computed  on  the  smaller 
diameter  of  the  converging  section.  For  included  angles  greater  than 
30°,  it  is  recommended  that  the  procedure  for  abrupt  area  decrease 
be  applied  to  converging  sections. 

3. 3  Bends 

Loss  data  for  bends  are  shown  in  design  chart  Dr  12,  in  the  form  of  a 
low  speed  (incompressible)  flow  loss  coefficient,  and  a  Mach  Number 


correction  for  high  speed  flow.  The  low  speed  loss  coefficients  have 
been  adapted  from  the  correlation  of  Henry,  Reference  8.  Mach  Num¬ 
ber  corrections  have  been  computed  from  the  data  Of  Young,  Reference 
16,  and  Higginbotham,  Reference  17. 

Few  investigators  have  determined  bend  losses  experimentally  for 
high  speed  subsonic  flows.  Some  general  trends,  hovvever,  are  ap¬ 
parent  from  the  results  available.  Data  of  References  16  and  17  show 
that  bends  of  circular  cross-section  tend  to  choke  at  an  inlet  Mach 
Number  of  about  0.7-0.75.  Near  this  critical  Mach  Number,  loss 
coefficients  increase  rapidly.  Bends  of  square  cross-section  tend  to 
choke  at  a  lower  inlet  Mach  Number,  about  0.  6-0.  7,  and  exhibit  simi¬ 
lar  increases  in  loss  coefficients.  Loss  coefficients  for  bends  depend 
to  some  extent  on  inlet  vei  u  it  ..istribution  (boundary  layer_thickness), 
but  this  effect  is  apparent!  .'  sn  iter  than  that  of  compressibility,  con¬ 
trary  to  the  results  for  ciru^  r» 

No  data  has  been  found  in  n:  i:„-raiure  for  Mach  Number  effects  in 
high  speed  flows  in  constant  area  bends  with  turning  vanes.  Reference 
18  reports  the  results  of  a  study  of  a  vaned  bend  of  rectangular  cross- 
section  with  an  area  increase  of  approximately  50%  (diffusing  bend), 
for  Mach  Numbers  from  0. 15  to  0.  40.  These  results  show  that  the 
loss  coefficient  for  such  a  bend  is  as  low  as  15%  of  the  loss  in  a  con¬ 
stant  area  bend  without  vanes,  with  the  loss  coefficient  increasing 
with  Mach  Number,  about  30  to  40%  greater  at  M  =  0.  5  than  the  in¬ 
compressible  loss.  This  percentage  increase  in  loss  coefficient  with 
Mach  Number  is  approximately  the  same  as  the  percentage  increase 
for  a  bend  without  vanes.  This  data,  and  the  data  of  Henry,  Reference  8, 
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show  that  large  reductions  in  bend  loss  coefficients  are  possible  by 
employing  vanes.  The  vanes  to  which  diese  data  apply,  however, 
are  airfoil  sections  (see  respective  References  for  design  information), 
and  apparently  require  considerable  sophistication  in  design  technique, 
to  achieve  low  losses.  For  such  vanes,  the  per  cent  increase  in  loss 
coefficient  With  Mach  Number  is  apparently  of  the  same  or-er  as  for 
bends  without  vanes. 


For  bends  with  simple  "splitters"" (i.  e. ,  vanes  which  are  not  of  airfoil 
^shape),  it  is  recommended  that  the  incompressible  loss  coefficients 
from  BuShips  Design  Data  Sheet,  Reference  1,  be  applied  in  conjunc¬ 
tion  with  the  Mach  Number  correction  curves  for  bends  without  vanes. 
There  is  no  specific  data  to  support  this  procedure,  except  the  indi¬ 
cations  from  Reference  18  that  percentage  changes  in  loss  coefficient 
with  Mach  Number  are  of  the  same  order,  regardless  of  whether  the  bend 
has  vanes  or  not. 


The  proper  method  of  computing  ducting  components  which  are  close- 
coupled  is  suggested  by  some  observations  reported  by  Young  and 
Green  in  Reference  16.  *"  Close  -coupled"  components  refers  to  two 
adjacent  ducting  components,  such  as  a  bend  followed  by  a  diffuser, 
or  a  90°  bend  followed  by  another  bend,  with  less  than  4  diameters  of 
straight  duct  separating  theni.  For  such  component  arrangements,  it 
is  recommended  in  the  calculation  procedures  that  follow  that  the  inlet 
Mach  Number  to  the  second  component  be  considered  the  same  as  the 
inlet  Mach  Number  to  the  first  component.  This  recommendation  is 
based  on  the  following  observations  of  Reference  16. 

Young  and  Green  have  reported  from  visual  observations  (Schlieren 
photographs)  of  high  speed  bend  flows  that  in  the  exit  piano  of  tlto 
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there  is  an  isentropic  core  flow  in  which  the  Mach  Number  is  unaffected 
by  losses  in  total  pressure.  In  this  exit  plane,  losses  are  concentrated 
in  the  low-energy  fluid  in  a  thickened  boundary  layer  near  the  bend  walls. 
This  total  pressure  loss  does  not  affect  the  Mach  Number  in  the  core 
flow  until  a  distance  of  about  4  diameters  downstream,  when  mixing  of 
the  core  flow  ar.d  the  boundary-layer  flow  has  been  completed.  From 
this,  it  appears  that  if  one  must  choose  a  one- dimensional  Mach  Number 
to  characterize  the  flow  in  the  exit-plane  of  the  bend  (as  is  necessary 
for  the  calculation  procedures),  the  proper  choice  is  to  assume  the  Mach 
Number  unchanged  from  bend  inlet  plane  to  bend  exit  plane.  Four  di¬ 
ameters  downstream  of  the  bend  exit,  after  mixing  has  occurred,  the 
effect  of  the  total  pressure  losses  in  the  bend  will  be  felt  all  across  the 
duct  and  the  one-dimensional  Mach  Number  will  have  increased. 

3.  4  Screens  and  Gratings 

The  subsonic  compressible  flow  through  round-wire  screens  or  gratings 
(screens  of  sharp- edged  elements  of  any  geometry)  has  been  considered 
by  Cornell,  Reference  12.  Cornell  gives  loss  coefficients  for  such 
screens  and  gratings  in  duct  flows,  for  various  screen  solidities  (ratio 
of  open  area  in  screen  to  ifs:  r.-um  duct  flow  area)  and  upstream  duct 
Mach  Number.  These  ha  .e  boon  taken  directly  from- Reference  12  and  arc 
included  as  design  chart  !'  13' 

The  loss  coefficients  of  C.  nn-.-'i  uo  not  apply  to  screens  at  the  mouth  of 
a  sharp-edged  (i.e.,  no  b  dimwUth)  intake  duct.  Cornell’s  calculations 
for  sharp-edged  screens  have  been  extended  for  screens  at  the  inlet  of 
a  sharp-edged  entrance,  based  on  an  extrapolation  of  the  data  for  con¬ 
traction  coefficients  given  in  Figure  7  of  Reference  12.  Results  are 
shown  in  design  chart  D-8. 


4. 0  CALCULATION  PROCEDURES  AND  EXAMPLES 


4.  i  Interpretation  of  Procedure  Diagrams 

Calculation  procedures  for  all  ducting  components  are  described 
graphically  in  the  diagrams  P-i  to  P-10.  These  diagrams  are  to 
be  interpreted  according  to  the  following  examples,  taken  from 
procedure  P-2,  i. 

1.  Quantities  which  must  be  known  as  initial  data  for  the  calcu¬ 
lation  are  shown  in  circles,  e.  g. 


2.  Quantities  which  a  re  determined  by  the  calculation  are  shown 
in  boxes,  e.  g. 

Re 

3.  Directions  for  calculation  are  shown  on  linfes  connecting  ele¬ 
ments;  e.  g. ,  the  sequence 


indicates  that  the  design  chart  number  D-2  is  to  be  entered 
with  quantities  whic..  arc  known  as  initial  data,  m/D  and  T  , 
and  Re  is  to  be  read  from  the  chart.  The  sequence 


nMMpMUW 
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indicates  that  4f,  which  is  a  result  of  previous  calculation,  and 
L/D,  which  is  known  as  initial  data,  are  to  be  combined  by  cal¬ 
culation  to  give  4f(L/D).  Here,  the  equation  required  for  cal¬ 
culation  is  obvious;  i.e.,  4f(L/D)  =  4f  x  L/D.  In  cases  where 
the  equation  is  not  obvious,  it  is  included  in  the  box  containing 
the  result,  as  in  the  sequence 


The  result  of  the  indicated  calculation  (i.  e. ,  the  "output"  of  the 
box)  is  always  shown  on  the  right-hand  side  of  the  equation. 

4.  In  some  cases,  alternate  calculation  procedures  are  indicated, 
as  in  the  sequence  (taken  from  procedure  P-10). 


This  diagram  indicates  tin. :  the  upper  path  is  followed  if  L/D  <30 


and  the  lower  path  if  L  D  >  30. 


-■  2  Examples 

.  Examples  illustrating  the  use  wf  each  procedure  diagram  follow  the 
diagram.  Each  box  contains  the  numerical  result  of  using  the  design 
charts  or  of  calculation  as  prescribed  by  the  procedure  diagram. 
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4.  3  Closely  Coupled  Components 


Area  changes,  bends,  and  screens  are  assumed  to  be  followed  by  4 
diameters  of  straight  duct  in  which  the  flow  mixes  and  stabilizes. 
When  two  components  are  not  separated  by  4  duct  diameters,  the  cal¬ 
culation  procedure  is  modified  as  follows: 


1.  The  pressure  loss  of  the  first  component  is  determined  in  the 
usual  manner. 


2.  The  pressure  loss  of  the  second  component  is  determined  by 
assuming  the  same  entry  Mach  Number  as  employed  for  the 
first  component.  The  pressure  loss  for  the  combined  compon¬ 
ents  is  then 


Po2 

'P02 

p  . 

*-  -J  1st  component 

LpoiJ 

2nd  component 

LPolJ 

3.  The  exit  Mach  Number  for  the  combined  components  is  deter¬ 
mined  by  calculating  the  mass  flow  parameter  at  the  exit  of  the 
second  component,  , 


“TfC 

V 

"Af 

_po2  A2_ 

exit 

_poi  Ai_ 

/V 

entrance 

i 

> 

_ 

,  (e02/poi> 

total  total 

Exit  Mach  Number  is  determined  directly  from  this  value  of  mass 
flow  parameter  through  chart  D-la.  This  procedure  is  illustrated 
in  Section  5  for  a  screen  followed  by  a  bend. 


PROCEDURE  P-1 

MACH  NUMBER  CALCULATION 


Total  pressure  known 


Static  pressure  known 


PROCEDURE  P-1 


MACH  NUMBER  CALCULATION 


Total  pressure  known 


Static  pressure  known 


Known: 


Mass  Flow 
Total  Temperature 
Duct  Diameter 
Total  Pressure 
Static  Pressure 


'  m  =12.  8  lb/sec 

T  =540  R 
o 

D=1  ft 
pQ  =14.7  psi 
p  =  14. 2  psi 


PROCKFXJRE  P-2.  I 


STRAIGHT  IXJCT  (INTAKE) 
% 


DUCT 

mat'l 


CALC 


!*_  .Pc 


'b- 
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PROCEDURE  P-2.  1 


STRAIGHT  DUCT  (INTAKE) 


PROCEDURE  P-2.2 


STRAIGH  T  DUCT  (EXHAUST) 


CALC 


PROCEIXJRE  P-2.2 
STRAIGHT  IXJCT  (EXHAUST) 


I 


Hoi 
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PROCEIXJRE  P-3.1 

nuc  r  INLET  LENGTH  WITH 
BELLMOUTI I  ENTRANCE 


rROCIilXJ  Rli  l>  .V.  2 


IXJCI*  INLKT  UiNC.  Ill  Will!  SHARP  KOGlil) 

iiN  i  ranck,  wm i  or  wm iou  r  soki;i£n 
AT  KNTRAN'CK 


—  k 
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PKOCblXJ  Kli  I*  .1.  2 

IXJCT  INLIiT  LKNG  I'll  WITH  SHARP  HIX 
liNTRANCl-;.  WITH  OR  WITHOU  T  SCRIilS 
AT  KNTRANCK 


0-9 


Pom 


'•/ 


PROCEDURE  P  4. 1 

ABRUPT  AREA  INCREASE  (INTAKE) 


H 

__ 

0-9  ^ 

■ai 

PROCEDURE  P-4.2 

ABRUPT  AREA  INCREASE  (EXHAUST) 


PROCEDURE  P-4.  1 


PROCEDURE  P-4.2 

ABRUPT  AREA  INCREASE  (EXHAUST) 


! 


PROCEDURE  P-5. 1 


D-/o 


ABRUPT  AREA  DECREASE  (INTAKE) 


M, 

L 

D-JO  . 

/n 

Pas. 

I, 

_ L 

os 

v  . -ELXJRE  P-5.2 

■-A-’PT  AREA  DECREASE  (EXHAUST) 
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PROCEDURE  P-5,  i 


PROCEDURE  P-5.2 

ABRUPT  AREA  DECREASE  (EXHAUST) 


.4. 
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PROCEDURE  P-6. 1 
DIFFUSERS  (INTAKE) 


0 


PROCEDURE  P-6. 1 
DIFFUSERS  (INTAKE) 


CHOOSE  NC*S  VALUE  OF  A/* 

AMO  repeat  calculation 


PKOCIilXJKK  P-6.2 
IXFFUSERS  (liXI  IAUS'1) 


PIUKJKIXJKI-  l»-7.2 


ASSUME 


M, 


CHOOSE  NEW  VALUE 
Or  Mt  A  NO  REPEAT 
CALCULATION 


»V7T 
Po,  A 


ASSUME" 


PROCEDURE  PS.  | 

SCREENS  AND  GRATINGS  (IN  I  AKIK) 


PROCEDURE  P-S.2 

SCREENS  AND  GRA  TINGS  (EXHAUST) 
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PROCEIXJRE  P -*i.  1 


IXJCT  INLE'TS,  AL  TERNATE  PROCEIXJRE 
FOR  BELLMOUTI I  ENTRANCE, 

il  >  30 
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PROCEIXJRE  P-9. 1 


IXJCT  INLH’rS.  ALTERNATE  PROCEIXJRE 
FOR  BELLMOUTII  EN  TRANCE, 

—  >  30 
O 


PROCEDURE  P- 


PROCEIXJRE  P-9.2.2 


IXJCT  INLETS.  ALTERNATE  PROCEDURE 
FOR  BELLMOUTI1  ENTRANCE. 


ENTER  P-2.1  or 
P-2.2  AT  POINT  A 


J 

] 


PROCEIXJ 


ROCEIXJKE  !»•«*.  2. 1 


UC  r  INLETS,  ALTERNATE  PROCEIXJRE 
OR  BELLMOUTII  ENTRANCE, 


PROCEDURE  P-9.2.2 

DUCT  INLETS,  ALTERNATE  PROCEDURE 
FOR  BELLMOUTII  ENTRANCE, 


1 


PROCKIXFKIi  I’-H) 

IXJCT  INLliT,  ALViiRNA  Hi  PROCli  DURli 
FOR  SilARP-  KIXiliH  LN  I'RANCIi 


ENTER  P-3.2' 
AT  POINT 


PROUilXFKIi  P- 10 


IXJCT  INLliT,  AL  I  I-  KNA  I'l*  PROC’lilXJRLi 
FOR  SIIARl’-l-IXllin  liNTRANCIi 


CALC 


MACH  HUMBER 
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FIGURE  D-2 


rel.  a  ty  ve  Roughness 


\ 

! 


-68- 


^  •  •%  %  4< 


J  -  J 


FIGURE  D-3 


F/?/ct/oa/  facto*,  a f 


LET  LENGTH  I  .OSS  COEFFICIENT  4  f 
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5.  0  APPLICATION  OF  CALCULATION  PROCEDURES  TO  A 
DUCT  SYSTEM  ' 

The  duct  system  shown  in  Figure  1  was  chosen  to  illustrate  the 
application  of  the  calculation  procedures  and  to  provide  an  example 
comparison  of  pressure  and  power  losses  for  different  duct  diame¬ 
ters.  The  results  of  the  calculations  for  three  different  duct  diame¬ 
ters  of  this  duct  system  are  tabulated  in  Table  5. 1.  The  pressure 
losses  and  power  losses  are  plotted  in  Figure  2.  Note  the  severe 
increase  in  calculated  losses  as  the  duct  diameter  decreases. 

The  power  loss  penalty  associated  with  the  duct  pressure  losses  has 
been  calculated  from  a  linearized  analysis  described  in  Section  5.  2. 
The  pressure  losses  are  so  large,  however,  that  the  linearized 
analysis  should  be  regarded  as  giving  only  approximate  estimates 
of  the  power  penalties. 
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TABLE  5. 1 


CALCULATED  PRESSURE  AND  POWER  LOSSES 
FOR  VARIOUS  DUCT  DIAMETERS 


INTAKE  DUCT  (with  flow  rite  of  15  lb/sec) 


Duct  Diameter  (in. ) 

18 

12 

9 

Pressure  Loss  (pal) 

0.1 

0.77 

1.64 

Power  Loss  (HP) 

8.5 

65.5 

141 

EXHAUST  DUCT  (with  flow  rate  of  16  Ib/sec) 


Duct  Diameter  (in. ) 

18 

12 

9 

Pressure  Loss  (psl) 

0.5 

4.05 

9.7 

Power  Loss  (IF) 

42.7 

233 

'5.88 

TABLE  5. 1  (Coned) 

PRESSURE  RATIO  AND  ENTERING  MACH  NUMBER  FOR  EACH  SYST 


|  18"  DUCTS 

12"  DUCTS 

Components 

Entering 

Mach 

Number 

Tout 

Pressure 

Ratio 

Geometry 

Entering 

Mach 

Number 

!  Total 
!  Pressure 
Ratio 

r  ... 
i 

t 

Geon 

INTAKE 

« 

Entrance 

0. 102 

1.00 

0. 237 

1.00 

Screen,  S  «0. 8 

0. 998 

0. 989 

Bend,  90° 

0. 998 

R/D  *1.0 

0.996 

R.  D 

Straight  duct 

0.998 

L/D-20.4 

0. 240 

0. 980 

0. 996 

L/D 

R.D 

Bend.  90° 

1 

0.999 

R  ,  D  ■  1 . 33  ! 

0. 246 

Diffuser 

0.248 

0. 987 

• 

Aj,  / 

1 

• 

, 

» 

2  e  < 

Straight  duct 

0. 102 

1.00 

L/D-7:9 

1 

1 

TOTAL  INTAKE 
SYSTEM 

0.993 

0. 947 

i 

i 

EXHAUST 

i 

; 

: 

Exit 

0. 178 

0.400 

Straight  duct 

0. 177 

0.991 

L/D  *23. 3 

0.370 

0. 935 

L/D 

Bend,  90°  j 

0. 175 

R/D  *0. 75 

0.360 

0.974 

R  D 

Straight  duct  j 

0.174 

0.997 

L/D  *  12 

0.347 

0. 967 

L/D 

Bend.  180°  | 

l 

0. 173 

0. 998 

R/D- 2. 66 

0. 345 

0. 993 

R/D 

TOTAL  EXHAUST1 
SYSTEM  | 

1 

t 

_ ^1 

.  0. 875 

i 

TABLE  5. 1  (Corn'd) 

)  ENTERING  MACH  NUMBER  FOR  EACH  SYSTEM  COMPONENT 


1 

• 

i 

! 

12"  DUCTS 

9"  DUCTS 

Entering 

|  Total 

Entering 

Total 

Geometry 

Mach 

Number 

t  Pressure 
Rstio 

» 

Geometry 

Mach 

Number 

:  Pressure 
Ratio 

Geometry 

i 

i 

R/D  *1.0 

0.237 

r 

i 

i 

1.00 

0.989 

,  0.99O 

* 

R  D  ■  1. 5 

0.460 

r*  ' 

i 

i 

0.993 

i  0.949 

-  0.988 

R.D-2.0 

L/D-20.4 

0.240 

0.980 

0.996 

L/D  -32.5 

R  D  *2.0 

* 

0.507 

;  0.820 

L  D  -44.7 

R  D-l.33  , 

« 

0.246 

0. 700 

0. 976 

i 

R  D-2.67 

! 

0.248 

0. 987 

• 

A2'  Al“2,25 

0.750 

0.939 

A2/A  -4.0 

| 

L/D-7.9  | 

t 

, 

i 

! 

,  0.947 

1 

1 

» 

1 

2  e  -  io° 

i 

•i 

• 

i 

j  0.700 

! 

• 

2©  -11° 

• 

i 

0.400 

• 

: 

. 

0. 68S 

< 

i 

L/D  *23.  3  ! 

0.370 

|  0.935 

L/D -37.25 

0.490 

0.  810 

L  D  *51 

R/D  *0. 75 

0.360 

0.974 

RD-0.75 

* 

0.  470 

0.975 

R/D-l.O 

L/D  *  12 

0.347 

0. 967 

L/D  =20 

0.  429 

i  0. 933 

L  D  -  28 

R/D  *2. 66  j 
I 

|  0. 345 

0. 993 

R/D -4.0 

0.  425 

:  0. 991 

i 

1 

R/D  -  5.  33 

I 

! 


i  0. 875 

i 


i 


0.728 
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5. 1  Sample  Duct  Pressure  Logs  Calculation 

5.1.1  Intake  Duct 

Mass  flow  m 

Temperatur  s  T  ^ 

Pressure  Pol 

Duct  diameter  D 

A 

Mach  Number  ;  P- 1  a 

mT7  _  ls^o 

P0,A  "  (14. 7)  (144)  (.442) 


ic/sec 
540°  R 
14. 7  ptia 
9  In. 


*  0. 460 


Bellmouth  Entrance  ;  P-3. 1 

/ 

4rWD)r  *  ':47 

Po/Po),  -  142 

4TL/Dlqulv  *  °'03 

4TWD)2  *  147 ' 0-03 


M2  *  0.  465 

p<>/p<4  - K41 

P»2/Pol)en,  " 


1.44 


'  0. 373 


0.993 


•90- 


Screen,  solidity,  s  -0. 8;  P-8,  i  (Note  use  of  tame  entering  Much 
s  Number  for  closo- coupled  components) 


M1  -0.460 


s-0. 8  Cpo-0.40;  1  -  <P02/P0i>  "‘O51 


^  ... 


w. 


949 


screen 


Bend,  90  ;  P-7. 1  (Note  use  of  same  entering  Mach  Number  for  close- 

coupled  components) 


—  -  -=t-  *  20 

D 


_L5 

.75 


To  -  540 


1.5 


.75 


-  0.460 


2.0 


Re  *  2.  i  x  10 


Cp0  -  0.092 


*B 


1.07 


V1 


^  CP0  *  0. 0985 
P02 

1-— •  .  0.0125 
Pol 


Po2 
Pol 


0.9875 


bend 


-91* 


x  Pol.  .  mf7  .  0.373 

PcLA  po2  Po2A  0.931 


-  0.401 


M2  “  0.S07 

Straight  Duct  ;  P-2. 1 
D  «*  9  in.  »0. 75  ft 
L  -  33.5  ft  L/D  *44. 7 
Re  -  2.1  x  106  4j  a  0  0181 

c/D*  0.0005/.  75  «  0.00067 
4 T  L/D  *  0. 809 


.  - -  — I 

•92. 

Bend,  90°  ;  P-7. 1 

JL.il!  .  2.67 

D  0.75 

M3  -0.7 
Re  -  2. 1  x  106 

Cp0  -  0.074;  Kg*  1.29  (Note:  probably  chokes)  K^  -1.0 
K.  Kr  Cp^  *  0. 0955 

f  b  O 


-0.0237 

po3 


•93* 


■  4.0 


2e  ■  li.i0 


M4  -0.75 


CpQ  -0.112 

Kq  *  2.0  (Note:  no  data  *  may  choke) 


KDCpo* 


M,  -0.75 
4 


mlTo 

po4A4 


0.501 


mYfo’  „  po4  A' 


0.501 


po4A4 


Po5A5 


(0.939)  (4.0) 


0.1335 


Mg  -0.121 


0.224 


•>4 


*0f 


f.» 


Pressure  Ratio  for  Intake  Duct 

Pos  ['i'  ^  / 

■p— -  -  (0. 931)  (0. 820)  (0.  976)  (0.939)  -  0.700 


^  1.2  Exhaust  Duct 


m  *  16  Jb/sec 
T0  -  1500°R 
Pg  *  14.7  psia 
D  «  9  in. 

A  *  0.442  ft2 

Mach  Number  ;  p-ib 


mTfo\  n- 

~^k  °- 


662 


Mg  ■  0. 688 


Straight  Duct  ;  P-2.2 


L  ■  38. 25 


m  m  16 
D  *  0.75 


Tq  «  1500  °R 


Re  *  1.1  x  106 


21.3 


J_  .  0.00067 
D 

4f- 0.0182 


-k  -51 


95- 


47—  ■  0.928 
D 


I. 


M,  -  0.688;  (47— 2H1)  -  0.225; 

«  n 


47  ~max  j  -0.225  +  0.928  -  1. 
D  l4 


153 


1.102 


M.  -0.490 


o5 

?o4 


0.  810 


1.36 


Bend,  180°  ;  P-7.2 


—  -  1.0 

D 

M4  -  0.  490 
Re  «  1.1  x  10* 


Cpn  -  0. 177 
o 

Assume  Mg  -  0.  470 


Kg  *  1.086'  KgCpo  =0.192 


1- 


04 


o3 


0.025 


'  04 


63 


0.975 


i 


*^6 


ml/T 


Po3^ 


0.379 


ml/T 


0.  389 


M4  *  0. 487 


Straight  Duct  ;  P-2.2 


L  =  21  ft 


~  *21.3 


Re  *  1.1  x  10 


*  0.00067 
D 

47  *0.0182 

L  21 
D  *  0.75  ' 


28 


47  ii  *  0. 510 
D 


M3  *  0.  470 
4f  Lfnax 


1.35 


47  Lmax 


1.40 


Jo  73 


*  1.35  +  0.51  -  1.86 


M2  *0.429 


1.50 


^3 

Po2 


0.  933 
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Bend,  180  ;  P-7.2 


—  =  5. 33 


Re  «  1.10  x  10 


Assume  M^  ■  0. 425 


CP, 


Kg  *  1.02 

«  1.63 

WPo 


0.0748 


po2 

1-  ~  -  0.0094 
pol 


P02 

Pol 


0.991 


m1^ 

polA 


0. 351 


min 


po2A 


0.355 


M2  *  0.  430 


Pressure  Ratio  for  Exhaust  Duct 


o5 

oi 


*  0. 045 


(0.810)  (0.975)  (0.933)  (0.991)  =  0.728 
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Total  Pressure  in  Exhaust  Duct  ;  P-1  a 


M<  «  0. 688 


P  -Ac 
1  5 


0.  583 


m  * 


A  * 


p _ (16)(  1500) 

05  (0.  422)  (0. 583)  (144) 


17.  5  psia 


16  lbs/sec 
1500°R 
0.  442  ft2 


24. 0  psia 


00- 


5. 2  Effect  of  Pressure  Loss  on  Turbine  Performance 


Assume  tur’  ine  nozzles  choked 


(1)  Increase  in  turbine  exhaust  pressure,  p^: 


-  ICO  - 


Multiply  by 
8  W 


8P, 


mf  v-1 

1  +  ~m“  Cp  V  To3 
a 


k-l/k 


po3  po2 

Let  r  =  -  -  ■■■■■■■■  ,  design  pressure  ratio 

P  P4  Pol 


Linearizing, 

x 

AH? 


m. 


=  AW  = - 1  1  + 
rn  x 


k-1 


CpT  To3  7 


m a  *T  P  -  ~i  -p 


k-l/k 


V  P  4 


AH> 


m 


a 


-CT(^) 

1  P  4 


(2)  Decrease  in  compressor  inlet  pressure,  pQ^: 
Assume  compressor  pressure  ratio  does  not  change 


P„q  =  P~o  =  P, 


°3  K°2  P  01  ’  P03  rnp( 


p  *01 


8Wx  I  mf\  k-i  P4  l  p4 

-  +  i^nr  ”TcpTo3iT rr  t 


-1/k 


.  1 


8P, 


ol 


P  \  p  °1 


p02i 


Multiply  by  pQl, 


8  W 


®P, 


+  1  + 


ol 


m. 


rn 


k-1 


,  k-l/k 


*r  cp  t03  nr 


r  P  1 
p  *ol 


1 


- 101  - 


Linearizing, 


Aif.  =  AW  .♦  i  +  _L 

rn  x  l  nv 

a  V  a 


Cp  “T  To3  7 


l  \  k-l/k/Ap 


T  p 


0  /i 


Numerical  Example 


Turbine  conditions: 


SHP  «  1100  IT 


Fuel- Air  Ratio  -  0.014  lb  fuel/lb  air 
Efficiency  -  0.  70 

Turbine  Inlet  Temperature  -  2060°R 
Design  Pressure  Ratio  -  6.  5 


Assume  Cp  =  0. 24,  k  =  1.  4 


n  _  ,  f  n  k-1  ^  1 

CT  1  +  ~m~ )  ^T  Cp  IT  To3  r 

/  V  a/  p 


k-l/k 


I  • 286  77£ 

(1 . 014)  (0. 70)  (0. 24)  (0.  286)  (2060)  (— )  -44 £- 

0. 0  odO 


Inlet; 


=  83. 1  HP/lb/sec 


A  IP  =  m.  C 


a  T  p 


m  =15  lb/sec 


CT  =  83.1  IP /lb /sec 


AH’  =  -(15) (83.  3) (0.3) 
Exhaust: 


t M  =  -  ma  C_AR 

a  T  p 


ma  =  16  lb/sec 

CT  =  83. 1  HP/lb/sec 

M*  =  0. 425 
4 


14.7  -  10.3, 
14.7 

=  -  375  H> 


23. 4  psi 

21.  2  -  14.7 
14.7 


-  (16)  (83.  3)  (0  .  405)  =  -  5  8  8  HP 


NOMENCLATURE 


'  2 

A  Cross-sectional  arec.^ft^ 

CpQ  Total  pressure  loss  coefficient 

D  Diameter  of  duct,  ft 

D  Equivalent  diameter  of  non-circular  cross-section  duct, 

6 

. _ A _ _ 

4  x  wetted  perimeter 

Kn  Pressure  coefficient  compressibility  factor  for  bends 

D 

Pressure  coefficient  compressibility  factor  for  diffusers 
K.  Bend  angle  loss  factor 

<p 

L  Duct  length,  ft 

\ 

L  Duct  length  to  attain  M  =1.0,  in  a  constant  area  duct,  ft 

max  b 

M  Mach  Number 

9 

R  Gas  constant,  lb  force  -  ft/lb  mass.-  °R 

Re  Duct  Reynolds  Number  ,  p  V  D/p. 

Re  Inlet  length  Reynolds  Number,  pVx/p. 

T  Temperature,  °R 

V  Velocity,  ft/sec 

c  Local  velocity  -  f  sccnd,  ft/sec 

T  Friction  factor .  .a'  oveloped  flow 


-lOf 


NOMENCLATURE  (Cont1  d) 


Friction  factor,  inlc:  laminar  length 
f2  Friction  factor,  inlet  turbulent  length 
gQ  constant  =  32. 2  lb  mass- ft/lb  force-sec  J 

k  Ratio  of  specific  heats 

m  Mass  flow  rate,  lb  mass/sec 

p  Pressure,  psf 

s  Screen  solidity  =  open  area  in  screen/duct  area 

x  Distance  from  duct  entrance,  ft 


Greek 

«  Duct  wall  roughness,  ft 

2  e  Angle  between  diffuser  walls,  degrees 

ix  Viscosity,  lb  mass/sec-ft 

3 

p  Density,  lbmass/ft 

4>  Angle  of  duct  bend,  degrees 


Superscripts 


* 


State  at  M  =1. 0 


NOMENCLATURE  (Concl-  d) 


Subscripts 

o  Stagnation  or  total  state;  condition  of  fluid  brought 

rest  isentropically 

1  State  upstream  of  component 

2  .  State  downstream  of  component 

tr  Transition  point  from  laminar  to  turbulent  flow 

State,  four  diameters  downstream  of  state  1 


2' 


Special  Nomenclature  for  Turbine  Performance 
Calculation 

»  • 

Cp  specific  heat,  Btu  °R 

mf  k-1 

turbine  power  constant  -  1  +  C  -g- 

ma  ^ 

IP  net  power  delivered  by  turbine  unit,  horsepower 

W  net  shaft  work,  ft.  lbs/lb.  of  air 

mf  mass  flow  rate  of  fuel,  lbs/sec 

mg  mass  flow  rate  of  air,  lbs/sec 

Tp  ^  design  pressure  ratio 
compressor  efficiency 

7Jrp 


turbine  efficiency 
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